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By means of Raman spectroscopy of liquid microjets we have investigated the crystal¬ 
lization process of supercooled quantum liquid mixtures composed of parahydrogen 
(PH 2 ) diluted with small amounts of up to 5% of either neon or orthodeuterium 
( 0 D 2 ), and of 0 D 2 diluted with either Ne or pH 2 . We show that the introduction of 
Ne impurities affects the crystallization kinetics in both the pH 2 -Ne and oD 2 -Ne mix¬ 
tures in terms of a signihcant reduction of the crystal growth rate, similarly to what 
found in our previous work on supercooled PH 2 - 0 D 2 liquid mixtures [M. Kiihnel et ai, 
Phys. Rev. B 89, 180506(R) (2014)]. Our experimental results, in combination with 
path-integral simulations of the supercooled liquid mixtures, suggest in particular a 
correlation between the measured growth rates and the ratio of the effective parti¬ 
cle sizes originating from quantum delocalization effects. We further show that the 
crystalline structure of the mixture is also affected to a large extent by the presence 
of the Ne impurities, which likely initiate the freezing process through the formation 
of Ne crystallites. 
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I. INTRODUCTION 


The study of crystallization in supercooled liquids is of prime importance for understand¬ 
ing the fundamental mechanisms of crystal growth. It is relevant to diverse research areas 
such as microstructural control in alloys engineering or the formation of ice in the atmo¬ 
sphere. In addition, it is of increasing interest in the context of the glass transitioiP, as the 
maximum rates of crystal growth may themselves provide insights into the glass-forming 
ability of supercooled liquid^^Hl Among the systems that can be used as test beds for our 
current microscopic understanding of crystallization in supercooled liquids are binary mix¬ 
tures. The crystallization of binary mixtures differs significantly from that of pure fluids, 
displaying a rich phase behavior in dependence of particle size ratio and composition. For 
example, understanding why for some compositions certain metallic binary alloys invariably 
form crystal phases that compete with the glass formatioiP is an open and fundamental 
question for understanding crystallization and its interplay with the glass-forming ability. 

The dynamics and crystallization of supercooled liquid binary mixtures in dependence on 
composition and particle size disparity have been investigated by classical molecular dynam¬ 
ics simulations of systems interacting through a simple Lennard-Jones (LJ) pair potentiaP®. 
However, systematic experimental studies of such simple atomic and molecular systems have 
so far remained out of reach. Here, we present experimental results, obtained by employing 
the liquid microjet technique in combination with Raman light scatterin^^, on the crys¬ 
tallization process of supercooled binary liquid mixtures of either pH 2 or 0 D 2 diluted with 
small amounts of Ne impurities. 

The present work was partly motivated by our recent investigations of the crystallization 
kinetics of supercooled PH 2 - 0 D 2 liquid mixture^. The most striking observation in those 
studies was the strong dependence of the crystallization rate on composition: starting with 
a pure pH 2 system, the crystal growth was found to slow down considerably with increas¬ 
ing amount of 0 D 2 , reaching a maximum in the crystallization time for the (pH 2 )gg ( 002)20 
mixture. This behavior may appear surprising because the intermolecular interaction po¬ 
tential is isotope independent!^. However, quantum delocalization due to zero-point motion 
increases the radius of the pH 2 and 0 D 2 molecules to a different extent due to their different 
masses, resulting in an effective 0 D 2 to pH 2 particle size ratio smaller than unity. This 
feature led us to interpret the observed growth rates in terms of packing effects linked to the 
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interplay between composition and particle size ratio in these quantum binary mixture^. 

Here we have extended the above studies to supercooled pH 2 -Ne and oD 2 -Ne liquid mix¬ 
tures, addressing in particular relevant structural aspects of the crystallization process in 
the Ne mixtures as well as in PH 2 - 0 D 2 mixtures not taken into account in our recent workP'. 
In principle, a Ne atom can be viewed as an isotopic impurity when mixed to pH 2 or 0 D 2 
because its interatomic spherical potential is very similar to those of pH 2 and 0 D 2 . Indeed, 
the LJ potential parameters (potential well depth e and interaction length a) for the H 2 -H 2 
interaction are a = 2.96 A and e = 34.2 K, while for H 2 -Ne they are a = 2.8745 A and 
e = 35.49 K (Ref. 11411 . However, Ne is 10 times heavier than pH 2 and, as a consequence, 
the magnitude of quantum delocalization is much smaller, resulting in a signihcantly smaller 
effective size when compared to pH 2 or 0 D 2 . The pH 2 -Ne and oD 2 -Ne liquid mixtures in¬ 
vestigated here represent thus ideal systems to further explore the effect of composition and 
particle size ratio on the crystallization of supercooled binary liquid mixtures. 


II. EXPERIMENTAL DETAILS 

The supercooled liquid mixtures were produced by the microjet technique and probed by 
Raman scattering, as described in Refs [TT] and HSl A schematic view of the experimental 
setup is shown in Fig. The liquid jet, which propagates at a speed v ~ \J2Pq/ p {Pq 
is the source pressure and p is the liquid density), rapidly cools well below the melting 
temperature until it crystallizes, producing a continuous solid hlament several cm lon^^. A 
crucial feature of our approach is represented by the correspondence between the distance 
along the jet propagation direction, z, and time, t = z/v, allowing for a probe of the 
crystallization kinetics with sub-microsecond time resolution. 

The liquid mixtures were continuously injected into vacuum through a 5 /im-diameter 
glass capillary nozzle. The glass capillary, mounted on a microactuators stage allowing a 
displacement of the entire nozzle assembly, and thus of the probed volume (Fig. [^, along the 
X, y, and z directions with an accuracy of better than 1 /rm, was cooled by a continuous flow 
liquid helium cryostat, and its temperature Tq was actively stabilized within ±0.1 K. The 
vacuum chamber was evacuated by a 2000 1/s turbo molecular pump providing a background 
pressure below 3 x 10“^ mbar. The nuclear spin variants pH 2 and 0 D 2 were produced by 
continuous catalytic conversion at 17 and 22 K, respectively, from 99.9999% and 99.9% 
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FIG. 1. Schematic view of the experiment, drawn to scale. The liquid micro]et (light blue) 
flows along the z axis. The exciting laser beam (green) polarized parallel to either the z or y axes 
propagates along the x direction. Raman scattering is collected along the y direction. The depicted 
sampled volume is the projection, onto the scattering plane, of the spectrograph entrance slit and 
of the active area of the CCD detector, showing the space resolution. 
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purity natural H 2 and D 2 gases, respectively, resulting in 99.8% and 97.5% purity pH 2 and 
0 D 2 , respectively, the rest being represented by odd-J molecules. The room temperature 
gas streams of pH 2 , 0 D 2 , and 99.998% purity Ne were then mixed at the specihc ratios 
by two mass flow controllers, one for each component, working at a minimum flow rate of 
20 normal-ml/min. The equilibrium solubility of neon in hydrogen in the liquid phase is 
limited to about 5%, and to about twice that value for neon in deuterium, whereas a phase 
separation occurs at higher Ne contenlP^^I^. The experimental conditions for all the mixtures 
investigated here are reported in Table [Tj 

The liquid microjets were probed by recording Raman spectra of the vibrational Qi(0) 
and rotational 5*0(0) transitions of pH 2 and 0 D 2 as a function of z. The Ne impurities could 
not be detected as they are Raman scattering inactive. Raman scattering was excited by 
a 4 W plane-polarized cw Ar+ laser beam at A = 514.5 nm, with an intracavity etalon for 
single-mode operation, and focused down to a diameter of ~ 14 pm onto the hlament. The 
polarization plane of the exciting laser beam could be rotated by means of a A/2 plate at the 
exit of the laser cavity. During the measurements of the rotational spectra the exciting laser 
beam was polarized along the y axis (Fig. in order to minimize the stray light resulting 
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TABLE I. Relevant experimental parameters (see the text for details) and impurity content for 
all the systems studied in the present work. 


System 

Impurity (mol %) Pq (bar) Tq (K) 

V (m/s) 

PH2 

- 

0 

8.5 

16.0 

142 

0 D 2 

- 

0 

13.0 

19.6 

115 


Ne 

0.3 

9.4 

17.0 

147 

pH2-Ne 

Ne 

1 

11.2 

17.5 

153 


Ne 

2 

10.8 

18.0 

137 


Ne 

1 

16.0 

20.0 

125 

oD2-Ne 

Ne 

2 

17.0 

20.0 

123 


Ne 

5 

11.8 

20.0 

95 


0 D 2 

1 

11.5 

16.0 

164 


0 D 2 

3.2 

8.9 

16.0 

143 

PH 2 - 0 D 2 

0 D 2 

4.6 

9.3 

16.0 

144 


PH2 

2.4 

16.3 

20.0 

130 


pH2 

4.3 

15.3 

20.0 

127 


from the elastic scattering, whereas for the vibrational spectra it was polarized along the z 
axis to record the Raman scattering in full. 

The Raman signal, which is proportional to the molecular number density, was collected 
at 90° with respect to both the laser beam and hlament axis, and was focused onto the 
20 /im entrance slit of the spectrometer by an optical system with x9 magnihcation. A 
supernotch hlter was used to block part of the very intense Rayleigh scattering from the 
liquid hlament. The spectrometer, with one meter focal length, was equipped with a 2360 
groove/mm grating and a back-illuminated CCD detector with 20 x 20 /rm^ pixels cooled 
by liquid nitrogen to 153 K. Its spectral resolution is ~ 0.11 cm“^ at 4150 cm“^ (pH 2 
vibrational band) and ~ 0.24 cm“^ at 180 cm“^ ( 0 D 2 rotational band). As illustrated in 
Fig.lU the space resolution along the jet axis was determined by the ~ 2 /rm projection 
of the spectrograph slit onto the microjet. Along the x radial direction 10 stripes, each 2 
pixels high, were read onto the CCD detector, allowing for a radial sampling of the jet. The 
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Raman spectra presented here were obtained by adding up just over the illuminated stripes. 


III. RESULTS 

A. Crystallization kinetics 

Examples of vibrational spectra are shown in Figs Ha) and l^b) for the ( 002)99 (Ne)j^ 
and (pH 2)99 (Ne)^ mixtures, respectively, evidencing the phase transition from the liquid to 
the solid. As mentioned above, the crystallization behavior of the Ne impurities could not 
be addressed during the present work because they are not Raman active. The vibrational 
spectra allow extracting the solidihed fraction of the sampled hlament volume. The evolu¬ 
tion of the 0 D 2 and pH 2 solid fractions measured for all the oD 2 -Ne and pH 2 -Ne mixtures 
investigated here are shown in Figsj^c) and|^d), respectively, where the axial distance has 
been converted into time as explained above. 

Overall, the crystallization kinetics features of Fig. are qualitatively similar to those 
reported recently for PH 2 - 0 D 2 mixture^, yet the crystallization slowdown observed upon 
addition of small amounts of Ne impurities is strongly enhanced. In Fig. id) we see that 
the presence of merely 0.3% Ne results in a signihcantly earlier start of the crystallization 
compared to pure pH 2 . In the case of the (pll 2)99 (Ne)^ and (pH 2 ) 9 g (Ne )2 mixtures the 
vibrational peak corresponding to solid pH 2 appears already within the hrst 100 pm distance 
from the orihce [Fig. |^b)], whereas the solidihcation process lasts about twice with respect 
to the pure pll 2 case. In the case of the oD 2 -Ne mixtures [Fig. |^c)] the presence of up to 
2% Ne does not actually lead to a clear earlier onset of freezing with respect to the pure 
case, though the subsequent crystal growth slows down with increasing Ne concentration in 
a similar fashion as for the pIl 2 -Ne mixtures. We recall that an earlier onset of freezing has 
been also observed in the PII 2 - 0 D 2 mixture^, especially upon addition of 0 D 2 impurities. 

B. Crystal structure 

The vibrational Raman spectra provide insights into the crystallization kinetics but they 
do not offer any straightforward information on the structure of the growing crystal. Such 
information can be retrieved from the rotational spectra, which we have probed for the 
pH 2 -Ne and oD 2 -Ne mixtures, as well as for diluted PII 2 - 0 D 2 mixtures. Since the rotational 
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FIG. 2. (a), (b) Selected normalized vibrational Raman spectra of 0 D 2 (a) and pH 2 (b) measured 
as a function of the distance z from the orifice (right scale) for the 1% Ne mixtures. The double¬ 
line shape of the 0 D 2 bands visible in (a) results from the presence of less than 3% of J = 1 pD 2 
molecules, with a 50-fold enhancement in the Raman scattering intensity with respect to the J = 0 
moleculeJ^. Panels (c) and (d) show the time evolution of the solid fractions extracted from the 
vibrational bands for 0 D 2 and pH 2 , respectively, in mixtures with mole percentages indicated on 
the right axis. The solid fractions range from 0 to 1, as indicated by the dashed lines on the left 
and on the right of the experimental curves, respectively. 

spectrum of the liquid is a broad band in all cases, we focus here on the structural properties 
after completion of the freezing process. The evolution of the rotational excitations during 
the liquid-to-solid phase transition in supercooled pH 2 has been discussed elsewher^. 

In Fig. I^a) and|^b) we present the rotational spectra measured for the PH 2 - 0 D 2 mixtures 
with 0 D 2 and pH 2 as impurity, respectively; the spectra were recorded at a distance of about 
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FIG. 3. Normalized rotational Raman spectra of 0 D 2 (blue) and pH 2 (red) in completely solidified 
filaments of diluted PH 2 - 0 D 2 mixtures. The impurity component is 0 D 2 and pH 2 in (a) and (b), 
respectively. The 0 D 2 and pH 2 impurity percentages are indicated on the right and left scales in 
(a) and (b), respectively. Vertical red and blue dashed lines indicate the position of the Raman 
lines measured for equilibrium hep crystals of pure pH 2 and 0 D 2 , respectively. The small peaks 
visible around 173 cm“^ in (a) and around 358 cm“^ in (b) are plasma lines from the Ar"*" laser. 


4 mm downstream from the orifice, corresponding to a propagation time of 40 /rs, long enongh 
for a complete crystallization in all cases (see Fig. [^. The rotational bands of both pnre 
solid pH 2 [Fig. [^a)] and 0 D 2 [Fig. [^b)] hlaments are characterized by fonr distinct peaks. 
This featnre has been attribnted to a random hexagonal closed packed (rhep) strnctnre, 
i.e., an alternating stacking of hexagonal closed-packed (hep) and face-centered cnbic (fee) 
crystal domain^. In fact, bnlk solid pH 2 and 0 D 2 slowly grown from the melt always 
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exhibit an hep structure, whose rotational Raman spectrum consists of three lines [the red 
and blue dashed lines in Fig. [^a) and[^b), respectively]. However, while the hep lattice is 
the most stable crystal structure, it was experimentally found that bulk solid pH 2 and 0 D 2 
grown fast from the vapor at very low temperatures also form the fee crystal lattice, which 
differs in energy from the hep lattice only by a factor of 10 ® (Ref. HSj). The rotational 
Raman band of the pH 2 fee structure consists of two peaks at ~ 350 cm“^ and ~ 356 cm“^, 
the latter overlapping with one peak corresponding to the hep structure. The addition of 
small amounts of the impurity leads to a broadening of the rotational features, as shown 
in Fig. I^a) and|^b), while maintaining the coexistence of the hep and fee lattices. This 
broadening of the rotational lines upon dilution with either 0 D 2 or pH 2 has been observed 
and interpreted previously as the result of the coupling between the rotational states in the 
solid mixture^. 

Turning now to the analysis of the rotational excitations relative to the impurity some 
striking features appear. For instance, we see in Fig. [^a) that at the lowest concentration 
of 1% 0 D 2 the 0 D 2 rotational spectrum exhibits a single line at ~ 179 cm“^. This line is 
very sharp and its wavenumber and width match those of the gas phase. This suggests that 
most of the 0 D 2 impurities occupy sites with no nearest neighbor 0 D 2 molecules, being able 
to rotate freely as a results of the larger average distance between the pH 2 molecules in the 
crystal matrix. With increasing 0 D 2 impurity concentration up to about 5%, we observe 
the gradual appearance of three additional sharp lines at the wavenumbers of the 0 D 2 bulk 
hep crystal latticd^. This can be rationalized in terms of nearest neighbor ( 002)2 pairs 
surrounded by pH 2 molecules, where the magnetic degeneracy of the rotational states is 
split along the intermolecular axis, giving rise to the hep triplelP^^. 

In the case of PH 2 - 0 D 2 mixtures with pH 2 as impurity [Fig. [^b)] we find a similar trend, 
yet with noticeable differences. The pH 2 rotational spectrum at the lowest concentration 
of 2% pH 2 is also characterized by a main peak whose wavenumber agrees with that of the 
free molecule, indicating that most of the pH 2 molecules are surrounded by 0 D 2 nearest 
neighbors. However, this Lorentz-shaped peak is significantly broader than in the case of 
isolated 0 D 2 impurities [Fig. [^a)] because the broader zero-point wave function of a pH 2 
molecule explore part of the repulsive potential of surrounding 0 D 2 molecules, affecting its 
rotational motion and thus giving rise to a broader Raman peak. The small side peak at 

352 cm“^, which grows with increasing pH 2 concentration, matches the lowest-energy 
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FIG. 4. Normalized rotational Raman spectra of 0 D 2 (a) and pH 2 (b) in completely solidified fila¬ 
ments of diluted pH 2 -Ne and oD 2 -Ne mixtures, respectively. The Ne impurity content is indicated 
on the right scale. 


excitation line of the hep triplet in the bulk (red dashed lines). This suggests again that a 
fraction of the pH 2 molecules in the lattice is arranged in nearest neighbor (pH 2)2 pairs, the 
remaining two hep peaks being probably hidden by the broad central peak. 

The rotational spectra of the solidified filaments containing Ne as impurity are presented 
in Fig. 1^ evidencing remarkable differences with respect to the case of the PH 2 - 0 D 2 mixtures. 
In Fig. I^b) we can see that the intensities of the hep lines of pH 2 at ~ 353 and ~ 356 cm“^, 
for a Ne concentration as low as 0.3%, drop to approximately half of their value for pure 
pH 2 , and become vanishingly small for the 2% Ne mixture. Simultaneously, the intensity of 
the fee line at 351 cm“^ increases, becoming the dominant feature, together with the line 
at 359 cm“^, at the highest Ne concentration of 2%. A similar trend is observed for the 
oD 2 -Ne mixtures [Fig. |^a)], though the two dominant peaks at ~ 175 and ~ 182 cm“^ in 
the 0 D 2 rotational spectra at the highest Ne concentrations are broader than those in the 
pH 2 spectra of Fig. |^b). Overall, the data of Fig. [^indicate that by adding small amounts 
of Ne to either pH 2 or 0 D 2 the crystal structure changes from an rhep configuration to a 
dominant fee structure, which is the equilibrium conhguration of solid neon. We point out 
that the slow crystallization from equilibrium H 2 -Ne and oD 2 -Ne liquid mixtures results in 
a physical separation into pure Ne and H 2 (D 2 ) solid phases, where the amount of the fee 
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domains is negligibly small at such low Ne concentration^^^®. 


IV. DISCUSSION 


At the microscopic scale the process of solidihcation from the melt is governed by crystal 
nucleation and the subsequent crystal growth. The overall picture that emerges from our 
crystallization kinetics data is that for the pure pH2 and 0D2 systems the limiting factor for 
crystallization is the formation of a critical crystallite; once a critical crystallite has formed, 
the subsequent crystal growth is fast (the growth rate for pure pH2 is about 30 cms“^). 
However, by adding small amounts of Ne atoms the nucleation time becomes shorter and 
the limiting factor for crystallization in this case is the reduced crystal growth rate. We 
present in Fig. a summary of the crystallization kinetics results by plotting the total 
crystallization time of the solvent, relative to the pure case, in dependence of the impurity 
content for the pH2-Ne and oD2-Ne mixtures, as well for the PH2-0D2 mixtured^. One can 
clearly distinguish two limiting cases, represented, on one side, by the barely signihcant 
increase of the crystallization rate of 0D2 upon dilution with pH2 (red squares) and, on the 
other side, by the strong dependence on the impurity content of the crystallization rate of 
pH2 upon dilution with Ne mixtures (green circles). 

According to the classical kinetic theorj^ the crystal growth rate is expressed in terms 
of thermodynamic and kinetic factors as 


u{T) 


k{T) 


_ g-AG(T)/fcsr 


( 1 ) 


where T is the temperature, AG{T) is is the difference in Gibbs free energy (per molecule) 
of the liquid and the crystal and is therefore the driving force for crystallization, k{T) is the 
crystal deposition rate at the liquid/crystal interface, and ks is Boltzmann’s constant. For 
the isotopic PH2-0D2 mixtures AG{T) can be computed from the experimental heat capacity 
data for the pure pH2 and 0D2 systems, and just a slight dependence on composition was 
founds. Due to the lack of experimental heat capacity data a similar direct determination 
of AG{T) was not possible for mixtures containing Ne. However, due to their nearly isotopic 
nature we do not expect that the thermodynamic factor plays a major role when describing 
the crystallization kinetics of pH2-Ne and oD2-Ne mixtures. Accordingly, the observed 
composition dependence of the crystal growth rate must be contained in the kinetic term 
k{T). 
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FIG. 5 . Crystallization time for microjets of diluted pH2-Ne and oD2-Ne liquid mixtures as 
determined from the pH2 (green circles) and 0D2 (green squares) solid fraction curves of Figsj^c) 
and [^d), respectively. Also shown are our previous data obtained for PH2-0D2 mixtured^ with 
0D2 (blue circles) and pH2 (red squares) as the impurity species. The plotted data indicate relative 
values with respect to the crystallization time measured for either pure pH2 or pure 0D2 liquid 
jets. The dashed lines are guides to the eye. 

The empirical Wilson-Frenkel model assumes that crystal growth is an activated process 
in which mass transport, expressed in k{T) through the diffusion constant, represents the 
most important limiting factor for the growth ratd^U Pow does mixing affect the trans¬ 
port properties, and consequently the crystallization kinetics, to the extent reported here? 
A number of simulation studies proposed that a lower particle diffusivity can be associated 
to dense packing effects resulting from the development in the liquid of some local coordina¬ 
tion geometric^. Recently, the composition dependence of transport properties observed in 
experiments with binary metallic liquids has been explicitly linked to corresponding changes 
in the density of packin^^^ES] indeed, from a pure geometrical point of view the packing 
efficiency in a mixture of hard spheres of specihc size ratio is a function of composition!^. 
In this respect, path Integral Monte Carlo (PIMC) simulations of supercooled PII2-0D2 liq¬ 
uid mixtures have evidenced that, as a result of the different effective sizes of the pH2 and 
0D2 isotopes, there is a slightly higher probability to hnd icosahedral-like order around an 
0D2 molecule than around a pll2 moleculd^. Among the possible local geometrical struc¬ 
tures icosahedra are those that are most densely packed and have been associated to the 
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FIG. 6 . PIMC simulations results for pH2-Ne and oD2-Ne with 3 % Ne computed at T = 13 K and 
T = 17 K, respectively, (a) Partial radial distribution functions g{r) representing the three pair 
correlations for the pH2-Ne mixture, (b) Probability distribution p{wq) for icosahedral-like order 
for the pH2 (red), 0D2 (blue), and Ne (green) particles. The green solid line refers to the pH2-Ne 
mixture whereas the dashed lines refers to the oD2-Ne mixture. The inset shows an enlarged view 
of the tail on a logarithmic scale, emphasizing the differences between the different species species 
at high negative values of the local bond order parameter w^. Note the similar probabilities for 
local icosahedral packing around a pH2 or 0D2 molecule. 


glass-forming ability of supercooled liquid^. 

We have performed analogous PIMC simulations of the Ne mixtures (see Appendix for 
details), and the main results are shown in Fig. In Fig. |^a) we plot the partial radial pair 
distribution functions computed for a (pll2)97 (Ne)3 mixture. A local bond order analysis (see 
Appendix) shows that the observed differences in the average static inter-particle correlations 
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translate into a higher degree of icosahedral-like order aronnd a Ne atom than aronnd a pH2 
molecnle, as shown in Fig. [^b). However, as a resnlt of the higher Ne mass compared to 
0D2 this tendency for local ordering aronnd a Ne impnrity in the pH2-Ne mixtnres is more 
prononnced than that fonnd aronnd an 0D2 molecule in PH2-0D2 mixtures. Furthermore, 
the comparison in Fig. |^b) of the distributions computed for the pH2-Ne and oD2-Ne 
mixtures shows also that the probability for local ordering around a Ne impurity in the 
former mixture is appreciably higher than in the latter mixture. These results thus further 
support a correlation between dense packing and slower crystallization kinetics, offering in 
addition a suggestive trend: the magnitude of the observed reduction in the crystal growth 
rate upon mixing appears to be roughly correlated with the effective particle size ratio for 
the specihc mixture. In fact, one can obtain an indirect estimation of the effective particle 
size ratio by comparing the hrst peak positions in the partial radial distribution functions 
computed along the PIMC simulations; it turns out that the ratio of the hrst peak positions 
in the pH2-Ne and pH2-pH2 radial distributions is about 0 . 966 , to which corresponds the 
strongest crystallization slow down (Fig. [^. This ratio becomes about 0.971 for the oD2-Ne 
system, whereas for the PH2-0D2 mixtur^^ the ratio of the hrst peak positions is about 
0.993 and 1.014 with 0D2 and pH2 as impurity, respectively. 

The rotational Raman spectra reported here oher additional insights into the crystal¬ 
lization process in supercooled quantum binary mixtures. We have seen (Fig. that the 
presence of small amounts of Ne impurities tends to force both pH2 and 0D2 to assume a 
crystal structure that can be associated to an fee lattice, which is the equilibrium conhgu- 
ration of solid Ne. It is unlikely that the presence of a tiny amount of Ne impurities at the 
1 % level and randomly distributed in the lattice would be able to modify the host pH2 or 
0D2 crystal structure to such an extent. The immediate interpretation of these observations 
is that the Ne impurities trigger the crystallization process with the formation of Ne-rich 
fee crystallites that act as nucleation sites and thus subsequently grow into the surrounding 
pH2- or oD2-rich supercooled liquid by keeping their initial crystal structure. The growth 
into an fee crystalline structure can be facilitated by the small energy difference between 
the fee and hep lattices. The formation of Ne-rich agglomerates is likely driven by compo¬ 
sitional fluctuations combined to the deep supercooling and the strong Ne-Ne correlation 
[Fig. [^a)]. That a clear tendency to assume a specihc crystal structure is not seen in the 
PH2-0D2 mixtures (Fig. is due to the fact that both pH2 and 0D2 crystallize into an hep 
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latticd^. This simple picture would also explain the earlier start of crystallization upon 
addition of the Ne impurities, especially evident for the pH2-Ne mixtures [Fig. [^d], as the 
probability for the formation of large impurity clusters increases with increasing Ne amount. 
The slight differences in the occurrence of the onset of freezing in the pH2-Ne and oD2-Ne 
mixtures [compare Figsj^c) and|^d)] may be partly related to the equilibrium D2-Ne phase 
diagranpSl^ which shows a slightly decrease of the liquidus temperature at small (< 3 %) Ne 
concentrations, in stark contrast to its increase observed in equilibrium H2-Ne mixtured^^^. 
As a final comment we note that a possible competition between the hep and fee structures, 
as suggested by our experimental data, might contribute to the slowdown of crystal growth 
in the pH2-Ne and oD2-Ne mixtures; the frustration of crystallization due to the conflict 
between different crystalline structures is an aspect that has been discussed in numerical 
studies of crystallization of binary mixtured^^ 

To conclude, we have shown that the combination of liquid microjets and Raman light 
scattering can provide highly valuable information on the crystallization process in super¬ 
cooled quantum liquid mixtures. Our results hint in particular at the importance of effects 
such atomic-level packing efficiency, related to effective particle size disparity originating 
from quantum delocalization, and crystal structure competition. It would be highly desir¬ 
able that the present results could be complemented by exploiting the capabilities offered 
by the liquid microjet technique in combination with state-of-the-art X-ray sources. Such 
studies would have the potential to shed further light on the physical origin of the findings 
reported here, and thus on the basic mechanisms of crystal growth. 
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Appendix: Simulation details 


We have simulated pH2-Ne and oD2-Ne mixtures with 3% Ne by assuming a simple 
LJ pair potential to model the interaction among the Ne atom^, the standard Silvera- 
Goldman potentiaP^ to model the interaction among hydrogen molecules, and the pair 
potential reported in Ref. |T 3 | to model the pH2-Ne and oD2-Ne interactions. We have 
simulated systems with N = 300 particles in boxes with periodic boundary conditions at a 
density given by the linear relation p = 0.97ppH2,oD2 + 0.03pNe, where ppHj = 0.0232 A“^, 
P0D2 = 0.0263 A“^, and pNe = 0.0373 are the densities close to triple-point conditions 
of the pure substances. 

Within PIMC the quantum system of N particles is mapped onto a classical system of 
N ring-polymer^. Each ring-polymer is composed of beads, whose number is fixed by 
expressing the full density matrix p{R,R',T) = (i?| exp{—H/kBT)\R') {H is the Hamilton 
operator, R and R' are many-body coordinates representing the positions of the N parti¬ 
cles, and T is the temperature), via a Trotter decomposition, as a convolution of higher 
temperature density matrices p{R, R',T'). In our simulations we adopted the pair-Suzuki 
approximation®! for p{R, R', T') with T' = 720 K and assumed a total number of 54 and 
41 convolutions for the pH2-Ne and oD2-Ne systems, respectively. The corresponding tem¬ 
peratures are T = 13.09 K and T = 17.14 K, respectively, which represent the estimated 
experimental average hlament temperatures of the pure system^^^, given the low Ne content, 
the same temperatures have been assumed also for the simulations of the mixtures. 

To determine the geometry of the focal environment (either crystalline or liquid) around 
a particular particle of the simulated quantum binary mixtures we used the identification 
schema based on the focal bond order (LBO) parameters introduced by Steinhardt et al.®^, 
modihed by averaging the LBO parameters over the nearest neighbors of the analyzed par¬ 
ticle and the particle itselP^. Since quantum delocalization invariably leads to large fluctu¬ 
ations in the beads’ positions, thereby deteriorating the LBO analysis, the LBO parameters 
were computed by using the centre of mass of each ring-polymer as the particle position®!. In 
addition, in order to guarantee a univocal choice of the nearest neighbors of a given particle 
we employed a three dimensional Delanuay triangulation. The averaged LBO parameters 
were then used to define different coarse-grained rotationally invariant parameters that are 
sensitive to different spatial symmetried®^. For example, the invariant Qg is ideally suited 
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FIG. 7 . Contour plots of the computed probability distributions p(Q4,Q6) for pH2 simulated in 
the bcc, hep and fee crystal phases phases at T ~ 12 K and p = 0.0240 A“^, and in a metastable 
liquid state at T ~ 13 K and p = 0.0232 A“^. Note that the use of periodic boundary conditions 
compatible with a specific crystal lattice allows simulating also crystalline states different from the 
equilibrium hep solid phase of pH2. 

to distinguish between the crystalline and disordered liquid-like configurations, whereas the 
invariant allows distinguishing between different types of crystal structures. Conversely, 
one can construct the invariant Wq that we used to analyze the tendency of the simulated 
systems to develop icosahedral-like structured. 

In order to prevent the crystallization of the mixtures during our PIMC simulations we 
adopted different strategies. We used a number of particles and a simulation box with 
side-ratios non compatible with close-packed crystal lattices; moreover, we generated the 
starting disordered conhguration for our simulations as reported in Ref. | 35 l The validity 
of our choices has been checked by computing the probability distributions p{Qi, Qe) for 
the ordered and disordered conhgurations, shown in Fig. for a pure pH2 systend, the 
probability distribution for disordered supercooled liquid pH2 is clearly well separated from 
the distributions computed for different crystal lattices, indicating the absence of any sign 
of crystallization in the conhgurations sampled for the metastable system. Similar results 
have been obtained for all the quantum binary mixtures considered in the present study. 
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